Generation and amplification of Raman Stokes and anti-Stokes waves 
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We present general analytical expressions of Stokes and anti-Stokes spectral photon-flux densities that 
are spontaneously generated by a single monochromatic pump wave propagating in a single-mode opti- 
cal fiber. We validate our results by comparing them with experimental data. Limiting cases of the general 
expressions corresponding to interesting physical situations are discussed. © 2008 Optical Society of America 
OCIS codes: 190.5650, 190.4380, 190.5890, 270.5530 



The third order nonhnearity seen by Ught propagat- 
ing in optical fibers is composed of two parts: a quasi 
instantaneous contribution coming from the electronic 
response, and a delayed Raman response coming from 
the coupling to molecular vibrations (phonons). The in- 
fluence of these two terms on the propagation of light are 
often considered separately, see for instance the treat- 
ment in [1]. But in many cases they must be treated 
collectively, together with the influence of dispersion 
which plays an essential role, particularly when the phase 
matching conditions are (even approximately) satisfied. 
That a holistic approach is necessary was first under- 
stood by Bloembergen and Shen [2], and confirmed ex- 
perimentally in later work. For instance it was shown 
that the gain of the Raman Stokes wave is strongly af- 
fected by the dispersion, and can even vanish when there 
is zero phase mismatch [3,4]. In addition it was shown 
in [2] , and confirmed experimentally in [5] , that the Ra- 
man anti-Stokes wave, which in the naive approach is 
exponentially damped, in fact grows exponentially at 
late times. This effect has important implications for su- 
percontinuum generation in photonic crystal fibers us- 
ing long pump pulses (ps and longer). In these systems 
spontaneous amplification of vacuum fluctuations at the 
anti-Stokes wavelengths can give rise to a blue detuned 
supercontinuum [6], see [7,8] for further studies and [9] 
for a review of supercontinuum generation. 

A quantum description of light propagation in optical 
fibers, including both the instantaneous Kerr non linear- 
ity and the Raman scattering has been developed [10] 
and used to quantify noise limits on squeezing [11], noise 
in x'-^'' parametric amplifiers [12], noise in coherent anti- 
Stokes Raman scattering [13], and noise in photon pair 
generation experiments [14, 15]. In the present work we 
apply this quantum theory, using the formulation of [12], 
to describe the spontaneous growth of Raman Stokes and 
anti-Stokes waves from vacuum fluctuations. Our ana- 
lytical predictions are compared to experimental results 
reported in [5,6] which had previously been fltted with 
semi-empirical formulae, or numerical calculations. 



We consider a continuous pump with complex ampli- 
tude ^/P e"^(^'*), where P is the power flowing through 
the flber and (j){x,t) = —ujpt + [kp + jPjx. In this ex- 
pression. Up is the pump frequency and kp is the (linear) 
wave number. Self-phase modulation gives a nonlinear 
contribution "fP to the total wave number. Here 7 is the 
third-order nonlinearity parameter of the fiber [1]. The 
quantum field operator 
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Aa{n,x) e"'"* I dfi + h.c 



describes the perturbations around this stationary solu- 
tion through the combined effect of Raman scattering 
and four- wave mixing. These are composed of symmet- 
rically detuned Stokes and anti-Stokes photons with re- 
spective frequencies ujs = ujp — n and Wa = ujp + f2, 
and wave numbers kg and ka- The operators As(Jl,,x) 
and Aa{i^,x) are destruction operators for Stokes and 
anti-Stokes photons. Their equations of motion can be 
deduced from the quantum theory of light propagation 
and solved analytically (so long as the perturbation field 
is small compared to the pump) [12]. 

Once As(ri,x) and Aa{^,x) are known, any physical 
quantity can be computed. Here we are concerned by the 
mean spectral photon-flux densities 
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Using the approach of [12], we find 
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The photon fluxes depend on three basic parame- 
ters: the pump power P, the detuning Q and the 
fiber temperature T. The pump power and the non- 
hnear response of the fiber are combined as x(^) — 
7P [(1 — /r) + /rXr(^)]j where the nonhnearity is de- 
composed into an instantaneous (electronic) part and a 
retarded (molecular or Raman) one. fr ~ 0.18 is the 
fraction of the total nonlinearity due to the Raman scat- 
tering and Xr-(^^) is the normalized spectral Raman re- 
sponse (Re[xr(0)] = 1, Im[xr] < 0). The complex param- 
eter K(n) [- {Ak{n)/2 f-Ak{n)xi^)]^/^ controls the 
growth rate of the Stokes and anti-Stokes waves. It de- 
pends both on x(^) on the linear phase-mismatch 
Ak{n) = ks + ka — 2kp. When the fiber is pumped suffi- 
ciently far away from the zero dispersion wavelength, the 
phase-mismatch Ak{n) is well approximated by /32^^^, 
where P2 is the group- velocity dispersion parameter. The 
real part of K{fl) is chosen positive so that it can be in- 
terpreted as the amplification gain for small signals at 
frequencies Wpztfi. Finally, the temperature T influences 
the photon fluxes through the factor 



n{n) = {exp{hn/kBT) - ly 



(5) 



related to the phonon population. 

The integral in Eq. Q can be carried out exactly, 
leading to somewhat cumbersome expressions. Here how- 
ever, we prefer to focus on two physically important 
limits : the spontaneous scattering limit {\k\x — > 0) 
and the stimulated amplification and wave-mixing limit 
(Re[K]x —f 00). The first one is of great importance for 
photon-pair generation, while the second one applies to 
supercontinuum generation. 

Considering the limit of small \k\x, Eqs. ^ give 
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up to the first order in \k\x. This contribution to the total 
photon flux is called the spontaneous Raman scattering. 
Stokes and anti-Stokes photons are emitted indepen- 
dently. A Stokes photon emission is always accompanied 
by the emission of a phonon of energy hCl, in contrast to 
the anti-Stokes process for which conservation of energy 
requires a phonon to be absorbed. For this reason, the 
anti-Stokes process in inhibited when the phonon pop- 
ulation is zero (low temperature limit). Using Eq. ([S]), 
one finds that ^iTn-\K.\x^o fa/ fs = exp[—hn/{kBT)] in 



accordance with the usual formulation theory of spon- 
taneous Raman scattering in fibers (see [5] and refer- 
ences within). Upon adding the extra term \x\'^ x'^ / {2tt) 
to Eqs. ([6K) and ([6b) one accounts for the spontaneous 
four-photon scattering process in which a Stokes photon 
and an anti-Stokes photon are emitted together after 
the absorbtion of two pump photons. This is the pro- 
cess used for photon-pair generation in fibers. Since it is 
only of second order in \k\x, spontaneous Raman scatte- 
ring always plays a detrimental role in photon-pair gen- 
eration experiments and is referred to "Raman noise" 
(see [14,15] for a more thorough discussion). It is inter- 
esting to note that linear dispersion has no impact on the 
values of the photon fluxes up to the second order in \k\x 
since they do not depend on Afc. However, the validity 
of the Maclaurin expansion is restricted to small values 
of Afc because |k| ^ (Afc/2) for fixed pump power and 
large Afc values. 

The asymptotic behavior for Re[K]a; ^ 00 is also sim- 
ple since we keep only the exponentially growing terms 
in Eqs. ^ : 
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Note that \K-{-i^ 



> \k — i^\ which implies that the 
Stokes wave is always stronger than the anti-Stokes wave, 
as expected. In Fig. [T] we have plotted the ratio fa/fs of 
the anti-Stokes to Stokes fluxes at the peak of the Ra- 
man gain {n/(2n) = 13.2 THz) at T = 300 K as a func- 
tion of -fP/Ak. In these circumstances, n — 0.14 and 
X = 7P (0.82 — i 0.25) for silica optical fibers. When 
7P/|Afc| ^ 1 pair creation dominates over Raman scat- 
tering in Eqs. (O, both in the normal and anomalous 




Fig. 1. Ratio of anti-Stokes to Stokes photon fluxes, at 
the peak of the Raman gain, as a function of jP/Ak. 
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Fig. 2. Evolution of the Stokes and anti-Stokcs spectral 
photon- flux densities with the pump power P. The pa- 
rameters are: 7 = 23 W^^km"^, a; = 2.85 m, /32 = 
50 ps^km"^ and T = 293 K. (Experimental data from 
[5].) 



dispersion regimes, and we have 

fs — fa 
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However fa/fs tends much faster to 1 in the anoma- 
lous dispersion regime than in the normal dispersion 
regime, as is apparent from Fig. [1] When the oppo- 
site limit 7P/|Afc| <C 1 is considered, the approximation 
K ~ |Im[x]| +i(Afc/2 4-Re[x]) +o{\x/Ak\) holds and 
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The ratio of anti-Stokes to Stokes fluxes is then given by 
fa/fs — IxP/A/c^. It is worth noting that it is indepen- 
dent of temperature. This simple ratio was first exhib- 
ited in [5]. However, it should be noted that its origin is 
more complicated than the argument given in [5] since 
the different processes of pair creation and phonon emis- 
sion/absorption all contribute to it. 

Even though Eqs. (O, ([8]) and ([9]) are valid for arbi- 
trary detunings Q, we now focus on ri/(27r) = 13.2 THz 
(corresponding to the Raman peak) in order to compare 
theory to experiments. 

Experimental Stokes and anti-Stokes spectral photon- 
flux densities in the 7P/|Afc| ^ 1 regime are available 
from [5]. They are plotted in Fig. [2] as a function of the 
pump power P. For the pump powers that are considered 
7P/|AA:| < 0.032. As seen from the figure, the asymp- 
totic formula ^ [dark grey curves] apply for P > 300 W 
only. At lower pump power the condition Re[K]a:: S> 1 is 
no longer satisfied. For very low pump powers, one can 
however apply the formula ^ [light grey curves] since 
|«;|a; <C 1. Outside these limiting cases Eqs. ^ reproduce 



very well the spectral photon-flux densities that are ob- 
served experimentally [black curves]. 

In [6], experimental data are obtained for a photonic 
crystal fiber with 7 = 150 W~^km~^ and P2 ^ 7 ps^/km 
(normal dispersion). Measurements have been carried 
out with a; = 3 m and P = 90 W, as well as with 
a; = 0.7 m and P = 400 W. In both cases Re[K]a:; > 1 
but since 7P/|Afc| is close to one neither Eqs. ([8|) nor 
Eqs. Q can be used to compute the ratio fa/fs- Our 
result ([7]) should be used instead. It predicts the ratio 
fa/fs = 0.028 for P = 90 W and fa/fs = 0.16 for 
P = 400 W (T = 300 K). These are in good agree- 
ment with the measured ratios [fa/fs = 0.016 and 0.22 
respectively) given the uncertainty on the experimental 
parameters in [6] and probable polarization effects. 

In summary we have developed an analytic theory to 
account for the growth of the Stokes and anti-Stokes 
waves from the combined effects of pair creation and Ra- 
man scattering. The results are in good agreement with 
earlier experimental observations. They should find ap- 
plications in the optimization of supercontinuum sources 
based on long pump pulses. 

This research was supported by the EU project QAP 
(contract number 015848) and the lAP programme, Bel- 
gium Science Policy, under grant P6-10. We thank S. 
Coen and D. A. Wardle for sharing with us their exper- 
imental data. 

References 

1. G. P. Agrawal, Nonlinear Fiber Optics, Third Edition 
(Academic Press, 2001) 

2. N. Bloembergen and Y. R. Shen, Phys. Rev. Lett. 12, 
504 (1964) 

3. E. Golovchenko, P. V. Mamyshev, A. N. Pilipetskii, and 
E. M. Dianov, IEEE J. Quantum Electron. 26, 1815 
(1990) 

4. F. Vanholsbeeck, Ph. Emplit, and S. Coen, Opt. Lett. 
28, 1960 (2003) 

5. S. Coen, D. A. Wardle, and J. D. Harvey, Phys. Rev. 
Lett. 89, 273901 (2002) 

6. S. Coen, A. H. L. Chan, R. Leonhardt, J. D. Harvey, J. 
C. Knight, W. J. Wadsworth, and Ph. St. J. Russell, J. 
Opt. Soc. Am. B 19, 753 (2002) 

7. A. B. Rulkov, M. Y. Vyatkin, S. V. Popov, J. R. Taylor, 
and V. P. Gapontsev, Opt. Express 13, 377 (2005) 

8. J. C. Travers, S. V. Popov, and J. R. Taylor, Opt. Lett. 
30, 3132 (2005) 

9. J. M. Dudley, G. Genty, and S. Coen, Rev. Mod. Phys. 
78, 1135 (2006) 

10. L. Boivin, F. X. Kartner, and H. A. Haus, Phys. Rev. 
Lett. 73, 240 (1994) 

11. J. H. Shapiro and L. Boivin, Opt. Lett. 20, 925 (1995) 

12. P. L. Voss and P. Kumar, Opt. Lett. 29, 445 (2004) 

13. D. Dimitropoulos, D. R. Solli, R. Claps, and B. Jalali, 
Opt. Express 14, 11418 (2006) 

14. Q. Lin, F. Yaman, and G. P. Agrawal, Opt. Lett. 31, 
1286 (2006) 

15. Q. Lin, F. Yaman, and G. P. Agrawal, Phys. Rev. A 
75, 023803 (2007) 



3 



